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ABSTRACT 

We report the detection of two very eccentric planets orbiting HD 4113 and HD 156846 with the CORALIE Echelle spectrograph 
mounted on the 1.2-m Euler Swiss telescope at La Silla. The first planet, HD 4113 b, has minimum mass of m sin? = 1.6 ± 0.2 Mjup, 
a period of P = 526.59 ± 0.21 days and an eccentricity of e = 0.903 ± 0.02. It orbits a metal rich G5V star at a = 1.28 AU 
which displays an additional radial velocity drift of 28 ms~Vyr observed during 8 years. The combination of the radial- velocity 
data and the non-detection of any main sequence stellar companion in our high contrast images taken at the VLT with NACO/SDI, 
characterizes the companion as a probable brown dwarf or as a faint white dwarf. The second planet, HD 156846 b, has minimum 
mass of rasin; = 1045 ± 0.05 Mjup, a period of P = 359.51 ± 0.09 days, an eccentricity of e = 0.847 ± 0.002 and is located at 
a = 1.0 AU from its parent star HD 156846 is a metal rich GO dwarf and is also the primary of a wide binary system (a > 250 AU, 
P > 4000 years). Its stellar companion, IDS 17 147- 19 14 B , is a M4 dwarf. The very high eccentricities of both planets can be 
explained by Kozai oscillations induced by the presence of a third object. 

Key words. Stars: planetary systems — Stars: binaries: visual — Techniques: radial velocities — Stars: individual: HD4113 — 
Stars: individual: HD 156846 



1. Introduction 

The CORALIE radial-velocity planet-search programme has been 
on-going for more than 9 years (since June 1998) at the 1.2- 
meter Swiss telescope located at La Silla Observatory, Chile. It 
is a volume-limited planet-search survey that contains aU main 
sequence stars from F8 down to KO within 50 pc and has a color- 
dependant distance limit for later type stars down to MO (Udry 
et al. 2000). Among the 1650 stars surveyed, 40 percent of them 
have a radial- velocity accuracy of 5 m s~^ or better and 90 per- 
cent of the sample is monitored with an accuracy better than 10 
ms"', the limitation being mainly due to photon noise. The re- 
maining 10 percent of the sample have a lower accuracy due to 
the lower signal to noise ratio and/or to the fast rotation of the 
targets. 

So far, CORALIE has allowed the detection (or has con- 
tributed to the detection) of 48 extra-solar planet candidates 
(e.g. Mayor et al. 2004; Segransan et al. 2007). This substantial 
contribution together with discoveries from various other pro- 
grammes have provided a sample of more than 230 exoplanets 
that now permits us to point out interesting statistical constraints 
for the planet formation and evolution scenarios (see e.g. Marcy 
et al. 2005; Udry et al. 2007; Udry & Santos 2007), and reference 
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therein for reviews on different aspects of the orbital-element 
distributions or primary star properties). 

In this paper we report the detection of two of the most ec- 
centric known planets , HD4113b and HD 156846b. Together 
with HD 80606 b (Naef et al. 2001) and HD 20782 b (Jones et al. 
2006), only four of the known planets have eccentricities larger 
than 0.8. The possible origin of these eccentricitie is still imder 
debate. Suggestions have been made for scenarios which allow 
the formation of eccentric orbits in a protoplanetary disk of plan- 
etesimals (Levison et al. 1998) or gas (Goldreich & Sari 2003). 
However, the presence of the third body in both HD 156846 b 
and HD4113b suggests that eccentricity pimiping might be at 
work in those two cases (Kozai 1962; Holman et al. 1997; 
Innanen et al. 1997; Mazeh et al. 1997). 

The paper is organized as follows. In Sect. 2 we briefly dis- 
cuss the primary star properties. Radial-velocity measurements 
and orbital solutions are presented in Sect. 3. In Sect. 4, we pro- 
vide some concluding remarks. 

2. Stellar characteristics 

2.1. HD4113, (HIP3391, SAO 192693) 

HD4113 is identified as a G5 dwarf in the Hipparcos cata- 
log (ESA 1997) and has an astrometric parallax, n = 22.70 ± 
0.99 mas which sets the star at a distance of 44.0 pc from the 
Sun. With an apparent magnitude V = 7.88 this implies an ab- 
solute magnitude of My = 4.66. According to the Hipparcos 
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Table 1. Observed and inferred stellar parameters for the stars hosting 
planets presented in this paper. Definitions and sources of the quoted 
values are given in the text. 



Parameters 




HD4113 


HD 156846 


Spectral Type 




G5V 


GOV 


V 




7.88 


6.506 


B-V 




0.716 


0.557 


n 


[mas] 


22.70 ±0.99 


20.41±0.93 


Mv 




4.660 


3.055 


7'eff 


[K] 


5688±26 


6138±36 




[cgs] 


4.40±0.05 


4.15±0.10 


[Fe/H] 


[dex] 


0.20±0.04 


0.22±0.05 


L 


[Le] 


1.22 


4.98 




[Me] 


0.99 


1.43 


V sin (' 


[kms^'] 


1.37 


4.45 


Age 


[Gyr] 


4.8-8.0 


2.1-2.7 



catalogue the color index for HD 41 13 is B-V - 0.716. Using a 
bolometric correction BC = -0.136 (Flower 1996) and the solar 
absolute magnitude Mhoi = 4.746 (Lejeune et al. 1998), we thus 
obtain a luminosity L = 1 .22 L©. Stellar parameters for HD 4113 
are summarized in Table 1 . 

A detailed spectroscopic analysis of HD 41 13 was performed 
using high signal to noise CORALIE spectra in order to obtain ac- 
curate atmospheric parameters (see Santos et al. 2005, for further 
details). This gave the following values: an effective temperature 
Teff - 5688 + 26 K, a surface gravity logg - 4.4 + 0.05, and a 
metallicity [Fe/H] = 0.20 + 0.04. Using these parameters and 
the Geneva stellar evolution models (Meynet & Maeder 2000), 
we derive a mass M = 0.99 Mq. According to the Bayesian 
age estimates of Pont & Eyer (2004), HD4113 is an old main- 
sequence star (4.8-8.0 Gyr). From the CORALIE spectra we de- 
rive V sin J = 1 .37 km s"' (Santos et al. 2002). 

2.2. HD 156846 (HIP84856 , HR6441, IDS 17147-1914 A) 

HD 156846 is a bright, metal-rich and slightly evolved GO dwarf. 
The astrometric parallax from the Hipparcos catalogue, n = 
20.41 ± 0.93 mas (ESA 1997) sets the star at a distance of 
49.0 pc from the Sun. With an apparent magnitude V - 6.50 
this implies an absolute magnitude of My = 3.055. According 
to the Hipparcos catalogue the color index for HD 156846 is 
B - V = 0.578. Using the method described in the previous 
section we derive a luminosity L - 4.98 Lq with a bolometric 
correction of BC = -0.049. Stellar parameters for HD 156846 
are summarized in Table 1 as well. 

A detailed spectroscopic analysis of HD 156846 gave the fol- 
lowing values: an effective temperature Tgff = 6138 ± 36K, a 
surface gravity logg = 4.15 ± 0.10, and a metallicity [Fe/H] = 
0.22 ± 0.05. Using these parameters and the Geneva stellar evo- 
lution models, we deduce a mass M = 1.43 Mq . According to 
the Bayesian age estimates of Pont & Eyer (2004), HD 156846 
is a moderately old main-sequence star (2.1-2.7 Gyr). From the 
CORALIE spectra we derive v sin j = 4.45 km s~' . 

HD 156846 is the primary star of a wide binary system 
(Washington Double Star catalog, hereafter WDS Worley & 
Douglass 1996) with an angular separaion of p = 5.1" and a 
position angle PA=75 deg. Based on its near infrared magnitude 
(J=9.405, H=8.92, ie. Cabrera-Lavers et al. 2006) and on the 
stellar evolutionnary models of the Lyon group (Baraffe et al. 
1998), its companion, IDS 17147-1914 B is identified as an early 
M dwarf of mass M - O.59M0. Using the positions given in the 
WDS, we obtain a projected binary separation of 250 AU. This 



translates into an estimated binary semimajor axis of 3 15 AU, us- 
ing the statistical relation a/r = 1.261 (Fischer & IVIarcy 1992). 



3. radial-velocity measurements and orbital 
solutions 

3.1. HD4113 

We took 130 spectra of HD4113 from October 1999 to October 

2007, yielding radial-velocity measurements with a typical 
signal-to-noise ratio of 50 (per pixel at 550 nm) with precision 
of ~ 3.5ms"'. 

Figure 1 shows the CORALIE radial velocities and the cor- 
responding best-fit Keplerian model with a linear drift. The re- 
sulting orbital parameters are P - 526.58 days, e = 0.90, 
K = 97.7 m s~\ imply the presence of a planet of minimum mass 
m sin J = 1.63 Mjup orbiting with a semi-major axis a = 1.28 AU. 
The orbital separation ranges from 0.12 AU at periastron to 
2.4 AU at apoastron. Residuals are somewhat larger than the pre- 
cision of the measurements {o-(o-c) - 8.4 m s '), and are proba- 
bly caused by a combination of stellar jitter and the inadequate 
modeling of the long period companion. 

In addition to the Keplerian signal, the radial-velocity mea- 
surements display a 28m s ' yr ' linear drift with a curvature 
of over -1.55 m s"' yr"^ over a time span of 8 years induced 
by a companion at a larger separation. Based on the duration 
and the amplitude of the drift, we are able to make some 
estimates on the nature of the outer companion of HD 4113. 
The lowest possible companion mass responsible for the drift 
has msinj = lOMjup which correspond to an orbital period 
of 11.5 year, an eccentricity of e = 0.5 and a radial-velocity 
semi-amplitude of 150 ms"'. However, this orbital solution 
is very unlikely since the outer planet-star separation would 
reach 2.5 AU at periastron which is the separation of the inner 
planet at apoastron. More realistic solutions therefore have 
longer periods. For instance, the shortest orbital period for a 
circular orbit is 22 years and correspond to a probable brown 
dwarf companion ( msinj = 14]V[jup) orbiting HD 4113 at 
a = 1.9 AU. The companion would reach the hydrogen burning 
limit at a separation of 20 AU with a period of 90 years. Radial 
velocities allow to conclude that the companion responsible 
for the observed drift is unlikely to be a planet (orbit stability 
reasons) but cannot discriminate between a brown dwarf and a 
stellar companion. 

To bring additional constraints on the outer companion's na- 
ture we obtained high contrast images of the target using the 
VLT and its adaptive optics system NACO in differential mode 
(Montagnier et al. 2007). We did not detect any object within 
a radius of 2.5 arcsec, which excludes any main sequence stel- 
lar companions to HD 4113 down to 0.2 arcsec ( about 8AU), 
(IVIontagnier et al., in prep.). We conclude that the companion re- 
sponsible for the observed drift is likely to be either a faint white 
dwarf like GJ 86 B (Mugrauer & Neuhauser 2005; Lagrange 
et al. 2006) that went undetected in our SDl images or it could 
be a brown dwarf located between 8 and 20 AU from the parent 
star with a period ranging from 20-90 years. 

3.2. HD 156846 

HD 156846 has been observed with CORALIE at La Silla 

Obervatory since May 2003. Altogether, 64 radial-velocity mea- 
surements with a typical signal-to-noise ratio of 75 (per pixel at 
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Fig. 1. Radial-velocity measurements as a function of Julian Date ob- 
tained with CORALIE for HD 41 13 superimposed on the best Keplerian 
planetary solution (top figure). The residuals are displayed at the bot- 
tom of the top figure and the phase folded radial-velocity measurements 
are displayed on the bottom diagram. 



550 nm) and a mean measurement uncertainty (including pho- 
ton noise and calibration errors) of 2.8 ms"' were gathered. 
Figure 2 shows the CORALIE radial velocities and the corre- 
sponding best-fit Keplerian model. The resulting orbital parame- 
ters are P = 359.50 days, e - 0.847, K - 464 ms"' , implying a 
minimum mass m2 sin / = 10.45 Mjup orbiting HD 156846 with a 
semimajor axis a - 0.99 AU. The orbital separation ranges from 
0.15 AU at periastron to 1.8 AU at apoastron. The orbital ele- 
ments for HD 156846 b are listed in Table 2. Residuals around 
the single planet orbital solution are larger than the precision of 
the measurements (cr(o_c) = 7.5 ms"') which points toward the 
stellar jitter since the the high eccentricity of HD 156846 b does 
not leave much room for a short period undetected planet. 

4. Discussion & conclusion 

In this paper, we have reported the detection of two very ec- 
centric planets orbiting two metal rich G dwarfs - HD 41 13 and 
HD 156846 - for which a third massive body is present in each 
planetary system. The origin of the wide range of eccentricities 



Fig. 2. Same as Fig. 1 for HD 156846. 



in extrasolar planets (represented in Fig. 4 as a e vs. logP di- 
agram) is still under debate, and different scenario have been 
proposed to explain the largest eccentricities such as the Kozai 
oscillations (Kozai 1962), chaotic evolution of planetary orbits 
in multiple systems and formation of eccentric orbits in a pro- 
toplanetary disk of planetesimals (Levison et al. 1998) or gas 
(Goldreich & Sari 2003). 

HD 156846 b and HD 4113 b are two of the four most ec- 
centric known planets, together with HD 80606 b (Naef et al. 
2001) and HD 20782 b (Jones et al. 2006) with eccentricities 
larger than 0.85. It should be noticed that the parent stars of 
these 4 planets are part of wide binary systems (Desidera & 
Barbieri 2007, for HD 20782 and respective discovery papers 
for HD80606, HD 41 13 and HD 156846) raising the question of 
the influence of the third body on the planets orbital parameters. 
It is therefore tempting to investigate the eccentricity pumping 
scenario a bit further by computing the different time-scales in- 
volved. 

The Kozai mechanism is effective at very long range, but its 
oscillations may be suppressed by other competing sources of or- 
bital perturbations, such as general relativity effects or perturba- 
tions resulting from the presence of an additional companion in 
the system. Regarding HD 4113 and HD 156846 , we have esti- 
mated their Kozai oscillation periods using Eqs (36) of Ford et al. 
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Table 2. CORALIE best Keplerian orbital solutions for HD 156846 
and HD 4113, as well as inferred planetary parameters. Confidence in- 
tervals are computed for a 68% confidence level after 10000 montecarlo 
iterations. Span is the time interval in days between the first and last 
measurements. cr(C* - C) is the weighted r.m.s. of the residuals around 
the derived solutions 



Parameters 




HD4113 


HD 156846 


y 


[km s"'] 


4.874 ± 0.05 


-68.540 ± 0.001 


slope 


[m s"' yr"'] 


27.8 ± 0.3 


- 


curvature 


\vr\ s"^ vr^^l 

Llll 3 J'l J 


-1.55 ± 0.17 




p 


[days] 


526.62 ± 0.3 


359.51 ±0.09 


K 


[m s-'] 


97.1 ±3.8 


464.3 ± 3.0 


e 




0.903 ± 0.005 


0.8472 ± 0.0016 


w 


[deg] 


-42.3 ± 1.9 


52.23 ±0.41 


To 


[JD-2.45 10«] 


3778.0 + 0.2 


3998.09 ± 0.05 


ai sin; 


[10-' AU] 


2.01 ± 0.05 


8.15 ±0.04 


f(m) 


[10-' Mo ] 


3.9 ±0.3 


558 ±9 


m2 sin ( 


[Mj„p] 


1.56 ±0.04 


10.45 ± 0.05 


a 


[AU] 


1.28 


0.9930 


N 




130 


64 


Span 


[years] 


8.0 


4.45 


O" (O-C) 


[ms-'] 


8.4 


7.7 



0.8 



0.2 
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Fig. 3. e VS. log P for the known extra-solar giant planets orbiting dwarf 
primaries. Filled dots are for planets orbiting single stars and open 
squares for planets orbiting a component of a multiple stellar system. 
Green stars correspond to stars hosting Netpune-like planets. HD 4113 b 
and HD 156846 b are located by filled red squares. 



(2000) . This yields to Pxazai = 5.5 10'' years for HD 156846 and 
Pkozoi = 6.2 10^ years for HD 4113 which indicates that both 
planet eccentricity can undergo Kozai oscillations with periods 
that are many order of magnitude shorter than the age of the 
systems. We also computed the apsidal precession period due to 
relativistic corrections to the newtonian equation of motion us- 
ing Eq. (4) of Holman et al. (1997). This yields Pcr = l.lGyr 
andf GR = 7.6Gyr for HD 156846 and for HD 41 13 respectively. 
This is also several order of magnitude larger that the Kozai os- 
cillation period. 

Provided that the orbital planes of the inner planet and the 
outer companion have different inclinations, the eccentricities of 
both HD4113 b and HD 156846 b could undergo Kozai oscilla- 
tions. That would be the most likely explanation for such high 
eccentricities. 
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